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ABSTRACT 

Using deep two-band imaging from the Hubble Space Telescope, we measure the color-magnitude 
relations (CMR) of E/SO galaxies in a set of 9 optically-selected clusters principally from the Red- 
Sequence Cluster Survey (RCS) at 0.9 < z < 1.23. We find that the mean scatter in the CMR in the 
observed frame of this set of clusters is 0.049±0.008, as compared to 0.031 ±0.007 in a similarly imaged 
and identically analyzed X-ray sample at similar redshifts. Single-burst stellar population models of 
the CMR scatter suggest that the E/S0 population in these RCS clusters truncated their star- formation 
at zi <z± 1.6, some 0.9 Gyrs later than their X-ray E/S0 counterparts which were truncated at zi ~ 2.1. 
The notion that this is a manifestation of the differing evolutionary states of the two populations of 
cluster galaxies is supported by comparison of the fraction of bulge-dominated galaxies found in the 
two samples which shows that optically-selected clusters contain a smaller fraction of E/S0 galaxies 
at the their cores. 
Subject headings: 



1. INTRODUCTION 

The color-magnitude relation (CMR) in galaxies in the 
cores of galaxy clusters is a fossil record of the integrated 
star-formatio n histories of the constituent galaxies (see 
lRenzinil[200l for a review). The slope of the CMR en- 
codes the build-up of the mass-metallicity relation, while 
its scatter is thought to reflect the v ariation in the forma- 
tion epochs of the cluster galaxie s dKodama &: Arimotol 
[j^lKauffmann fc CharlotHl99fit Bernardi et al.ll2005l ). 

Mounting observational data has culminated in a 
picture of galaxy formation that attempts to explain 
the conversion of star-forming "blue-cloud" galaxies 
into passively-ev o lving "red-sequ ence" galaxies (e.g., 
IFaber et al.ll2007t IMei et all 120091 ) . This picture cites 
some combination of the quenching of star formation in 
the blue-cloud and the merging of already red galaxies 
as driving the build-up of the red-sequence. 

The high resolution imaging delivered by HST/ ACS 
plays an important role in disentangling these effects, 
as it allows color measurements even for faint galaxies, 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope and obtained from the data archive at the Space 
Telescope Institute. STSci is operated by the Association of 
Universities for Research in Astronomy, Inc. under the NASA 
contract NAS 5-26555. The observations are associated with 
program 10496 

2 Department of Astronomy and Astrophysics, University of 
Chicago, Chicago IL 60637, USA 

3 Kavli Institute for Cosmological Physics, The University of 
Chicago, Chicago IL 60637, USA 

4 Visiting Associate, Observatories of the Carnegie Institution 
of Washington, Pasadena, CA 91101 

5 Department of Physics and Astronomy, University of Water- 
loo, Waterloo ON, Canada, N2L 3GI 

6 Department of Astronomy and Astrophysics, University of 
Toronto, Toronto, ON, Canada, M5S 3H4 

7 E.O. Lawrence Berkeley National Lab, 1 Cyclotron 
Rd.,Berkeley,CA,94720 

8 Department of Physics, University of California Berkc- 
ley,Berkeley,CA,94720 

9 Department of Physics and Astronomy, University of Utah, 
Salt Lake City, UT,84112 



and robust morphological measurements, which include 
either manual classifi cation of E, SO, Sp, and Irr (e.g., 
iPostman et al.ll2005ft or auto mated separation of E/SO 
types from later types (e.g.. [Abraham et al.l l200"7h . In 
particular, combined color and morphological data can 
be used to measure the color scatter of the E/SO galax- 
ies in the CMR, which places constraints o n the forma- 
tion e pochs of stars in cluster galaxies (e.g. jBower et all 
fl992T) . 

Until now, z ~ 1 systematic space-based studies of 
the CMR in the cores of galaxy clusters have been con- 
ducted almost exclusive ly using a handful of X- r ay se- 
lected galaxy cluster s (Ivan Dok kum etal] 1 19981 fe OQO; 
Blakeslee etaLI 12003b IMei et al.l 120061 : iBlakeslee et al.l 
20061: IMei et al]|2009| ) . The consensus picture from these 
studies is that stars in these galaxies were formed at 
z>2. 

In this Letter, we present the first precision statistical 
measurements of the scatter of the CMR in an optically- 
selected sample of galaxy clusters at z ~ 1. Where neces- 
sary, we adopt a flat, fl m — 0.30 cosmology with H = 70 
km s _1 Mpc -1 . 

2. DATA 

Imaging for this study was acquired as part of the 
HST Cluster SN Survey Program (Program Number 
10496, PI: Perlmutter). Nine optically-selected clusters 
were imaged at various telescope roll angles for a total 
of 5000- 16000s with the Advanced Camera for Surveys 
(ACS) on multiple visits (375s-500s per sub-exposure) 
in z SRn (F850LP) an d for 1000-8000s in i 775 (F775W) 
(|Dawson et al.l I2009D . Eight clusters are derived from 
the R ed-Sequence Cluster Survey fRCS. lGladders &: Yed 
12005ft . and the ninth, XLSS J0223.0+0436, was originally 
included in this samp le as an IR-detected cluster from 
the S pARCS survey (|Muzzin et alj l2008t IWilson et ail 
12008ft, but had previou sly been de tected in both th e 
IR (|Andreon et al.ll2005f ) and X-ray (jPierre et al.ll2004ft . 
We consider it here as an optical-IR detected cluster. 
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With the exception RCSJ J0337.8-2844 an extensive 
spectroscopic program has obtained 5-30 confirmed spec- 
troscopic members per cluster, e nough to compute a 
velocity dispersion in many cases (IGilbank et al.1 [20081 : 
iHicks et al.ll200l iGilbank et alJl2009f ). Three oflheX- 
ray se lected clusters, anal yzed most recently inlMei et aTl 
( 2009) are taken from (|Rosati et all 1 19981). and the 
WAR PS cluster was first presented in (jPerlman et all 
l200l . 

I mages for each cluster we re coadded using Multidriz- 
zle ()Koekemoer et all l2002f) . with a Gaussian kernel, 
pixf rac = 0.8, and a final output pixel size of 0.03 arc- 
sec/pixel. The detailed photometric and morphological 
pipelines de veloped for thi s work are described and tested 
in detail in iKoester et a l. ( 2009) . In brief, object detec- 
tion was performed with SExtractor (|Bertin fc Arnoutsl 
1996), and the resulting catalogs were manually cleaned 
of diffraction spikes, pupil ghosts, and saturated stars. 
The final gal axy catalog was d i stilled from the full cata- 
log following [Lcauthau d et al.l (j2007) . Total magnitudes 
are taken a s MAG-AUTO magnitudes with a correction ap- 
plied (e.g.jBernstein et~alll2002allb l; iBenitez et alj|2004: 
IBlakeslee et al.1 12006D appropri ate for E galaxies. As 
noted in iBlakeslee et al.l (|2003h . the differential smear- 
ing of the PSF between 1775 and 2:350 makes the PSF 
approximately 10% broader in z 85Q . We correct this us- 
ing a PSF matching scheme that places both the 2775 and 
Z850 images on the same footing by convolving the latter 
image with the PSF of the former, and vice-versa. 

The PSF model was built in each bandpass using 
the SExtractor MU_MAX - MAG_AUT0 parameters to select 
acceptable PSF stars in the coadded images, which 
were then visually checked. In a circularly symmet- 
ric ~ 1' radius, a single PSF model adequately de- 
scribes the aggregate distortion realized in the final im- 
age due to HST/ACS optics and tempor al variations 
thereof dJee et alj|2007t iRhodes et al.ll2007D . charge dif- 
fusion dKrist et al.l I2003D. the so-ca lled red-halo effect 
(|Sirianni et all 120051 : Uee et al.l I2007T ). and the drizzling 
process. For the precision photometry required in this 
study, this necessarily restricts us to considering the in- 
ner 500ft. -1 kpc of each cluster, where the PSF is ap- 
proximately constant, and the convolutions required for 
PSF-matching are minimally position dependent. 

In the coadded and convolved Z850 image of each clus- 
ter, we construct adaptiv ely-determined apertures with 
GALFIT ([Peng et al.ll2002D by fitting Sersic profiles con- 
strained to 1 < n < 4, thus returning an estimate of 
an effective radius, R e for each object. Nearby objects 
(< 3R e ) are simultaneously fit, and all remaining ob- 
je cts within 6R P are ma sked, and the sky estimated as 
in Haussle r et al.l (2007). The same aperture is applied 
in the 1775 band to measure the «775 — z%§o color. Er- 
rors on the aperture magnitudes are determined for each 
galaxy by comparing the repeat observations of objects 
as a function of magnitude. Typical errors range from 
0.02 mags at ZS50 — 21 to 0.05 magnitudes at zsso = 24. 

We use a n au tomated scheme, MORPHEUS 
([Abraham et al.1 l2007h to separate E/S0 galaxies 
from the irregulars and spirals, taking objects with 
a Gini coefficient > 0.45 and asymmetry parameter 
< 0.1 to be E/S0 galaxies. Because Zs50 stays redward 
of the 4000A break until z ~ 1.3, the signal-to-noise 



ratio (S/N) for early- type galaxies remains considerable. 
Based on comparison to the manual classifi cations for 
RDCS J1252. 9-2927 in iPostman et~aT1 pOOSl ). the E/S0 
galaxies we classify at z = 1.24 comprise a 90% pure 
and complete sample. 

3. MEASURING THE COLOR-MAGNITUDE RELATION 

Color-magnitude diagrams were created for each clus- 
ter, and the red-sequence color limits were determined 
by eye. Objects falling within these broad limits (~ 3cr 
of the final best-fit line) with luminosities greater than 
0.3L* (-Z850 = 24 at z = 1.24) were taken to be the red- 
sequence. The luminosity limit is strictly chosen to be 
bright enough to ensure the morphological classifications 
are robust. A line was fit to E/S0 objects within these 
bounds using a simple least-squares method, and was 
sigma-clipped at 3 a in i — z color. This was iterated un- 
til convergence to derive the slope and zero point of the 
best-fit line, and the scatter of E/S0 objects (including 
the rejected outliers) about thi s line was comput ed using 
a bi- weight estimator of scale (|Beers et al.lll990[) . 

The best-fit slopes of two clusters, RCS J2156. 7-0448 
and RCS J2345. 4-3632, are not robust to the color cuts 
used, due largely to the relative paucity of E/S0 galaxies 
at the faint end of their CMRs (z 8 50 > 22.5). We thus 
fix the slope for these clusters at —0.05 and the reported 
scatter is measured about this line; the resulting scat- 
ter about the best-fit line varies by only ~ 0.015 mags 
depending on the line used. 

To assess uncertainties, we run 1000 bootstrap resam- 
ples on each cluster, and compute the mean and standard 
error of the red sequence fits to these resamples. There is 
a further contribution to scatter due to photom etric er- 
ror. Using the method of IStanford et al.1 (|1998f ) we find 
that it contributes no more than 0.01-0.02 magnitudes 
in quadrature to the measured scatter; we correct the 
measured scatters accordingly. 

Finally, we note that the derived CMR is measured 
inside a fixed 0.5 h~ x Mpc cluster-centric radius, the ap- 
proximate observed field-of-view of most of these obser- 
vations. We do not attempt to compute -R200 or any 
other scale radius, which would require additional dy- 
nami cal data, X-ray d ata, or possibly wider field images 
(e.g.. lRozo etal"1l2008t L 

4. THE CMR IN z ~ 1 RCS CLUSTERS 

Table 1 shows the measured CMRs for our 9 optically- 
selected clusters and also our own measurements of a 
subset of archival X-ray selected clusters also measured 
in ^775 and z 850 . The Table also displays N E / S0 , the 
number of E/S0 galaxies used in measuring the CMR, 
and a richness -/V0.5, which is the number of red-sequence 
galaxies brighter than 0.3L*, within 0.5/1" 1 Mpc of the 
BCG. RCS J0337.8-2844 and RCS J2156.7-0448, the 
two least rich clusters, exhibit the largest scatter. RCS 
J2319.8+0038 and RCS J0439 6-2904, which are bo th X- 
ray detected (jHicks et al.ll2008l: IGilbank et al.ll2008h . and 
the X-ray selected clusters have relatively high A0.5 rich- 
nesses. 

In this Letter, we focus on the distribution of scatters. 
In the observed frame, it is apparent that the optical 
and X-ray samples have discordant average scatters in 
observed 1775— 2sso : the optical sample, at 0.049±0.008 is 
nearly 50% larger than the X-ray scatter at 0.031 ±0.007. 
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A i-test of the means of the two samples rejects the null 
hypothesis that the two samples come from the same 
underlying population at the 95% level. 

The scatter is though t to reflect variati ons in galaxy 
formation time (but see lMenci et al-lfeOOSh . It has thus 
been modeled using Bruz ual-Charlot single-burst models 
(jBruzual fc Ch ariot 20030 , which posit that all the stars 
in a given galaxy are formed in a coeval burst at some 
rcdshift, zt, and evolve passively thereafter. The scatter 
in the CMR results from the variation in formation times 
for the constituent galaxies that have formed in single 
bursts between reionization, z r , and some lower redshift, 
zi- 
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Fig. 1. — Bruzual-Charlot Single-Burst Models of the Scatter. 
Given the cluster redshift, z a i, 3 , stars with error bars represent 
the scatter of the red-sequence in X-ray selected clusters, while 
the normal crosses and errors show the scatters for the optically- 
selected clusters. The best-fit Bruzual-Charlot single-burst models 
for the X-ray and optical samples are shown as dashed and solid 
lines, respectively. Inset: \ 2 minimization for the X-ray and optical 
models, indicating the lowest redshift, zi, urat , of the epoch during 
which galaxies formed. The two samples clearly favor different 
minima in the context of this model. Alternatively, taking z^ urat = 
1.6 in the x-ray sample, the x 2 = 6.1, a result that is discordant 
with the optically-selected best-fit zi, urst = 2.1 model at the few 
sigma level. 



Using t he ACS ^775 and zsso response curves as de- 
scribed in iSirianni et al.l (|2005f ). we determine the ex- 
pected observed scatter in 1775 — z§5o using a single-burst 
spectrum, following the Monte Ca rlo method first de- 
scribed in Ivan Dokkum et alj (|1998h . Assuming a galaxy 
formed at some redshift, Zf, we create 10,000 "galaxies" 
by drawing from a random uniform formation age distri- 
bution corresponding to a redshift range of zi < Zf < z r . 
For each zi, the last redshift at which galaxies formed, 
the model predicts a scatter in 1775 — z 850 that would be 
observed at some redshift, z b s . For a given set of clus- 
ters, we then choose the best-fit model by minimizing 
X 2 : 



X 



= s 



T obs\2 



(1) 



where the "mod" and "obs" superscripts refer to the pre- 
dicted model and observed data scatters, and the 
error on each observed a ohs . In Figure 1, we show the re- 
sult of this process for both the optical (solid) and X-ray 



(dashed) samples. The optical sample favors zi = 1.6, 
and the X-ray sample favors z\ — 2.1, suggesting that 
the E/S0 galaxies on the red-sequence in X-ray selected 
clusters ended their star formation some 0.9 Gyr earlier. 

5. DISCUSSION 

The small number of E/S0 galaxies, N E / S0 , in some of 
the less rich clusters highlights a possible concern. Sim- 
ple Monte Carlo simulations of random draws from a 
normal distribution indicate that at iV = 10, even the 
bi-weight estimator of scale returns a scatter that is bi- 
ased high by 10%; this is reduced to 3% by N = 20, and 
< 1% at N=40. Contamination by non-cluster mem- 
bers also contributes to scatter. Randomly adding a sin- 
gle galaxy that passes the red-sequence color cuts adds 
as much as an additional 10% bias at N = 10. For- 
tunat ely, our analysis of the GOODS 9 (jGiavalisco et al.1 
l2004h fields indicate that the average background (for 
morphologically selected galaxies passing our color and 
luminosity cuts) is < 1 galaxy. 

Another concern is line-of-sight projection of large 
structures. RCS J0439. 6-2904 is known to show two red- 
shift peaks at z = 0-943 5 and z = 0.9681 (jGilbank et all 
I2007t ICain et al.l l2008f ). for a separation of ~ 7000 
km s _1 . A cursory ana lysis of the halo ca talogs from 
the VIRGO simulations (|Evrard et al.ll2002l ) shows that 
halo-halo projections of this sort are expected in about 
~ 10 — 20% of RCS clusters in this mass range at 
z ~ 1; a similar va lue is found in spectroscopic data 
(|Gilbank et al.ll2C)07h . At such a separation, Monte Carlo 
simulations show that two projected halos each with 
a v = 1000 km s^ 1 have a resultant measured scatter in- 
flated by only 0.003 magnitudes, which has no significant 
effect on our results. 

Spectroscopy does not yet exist for RCS J0337. 8-2844 
and spectroscopic confirmation is pending. Only lim- 
ited spectroscopy exists for RCS J2156. 7-0448 (5 con- 
firmed spectroscopic memb ers), and its meas ured X-ray 
flux is consistent with zero (Ifficks et al.ll2008ft. alt hough 
the strong lensing in this system (|Gladders et al.|[2003T) 
indicates that its mass is cluster scale. Regardless, we 
find that removing one or both of these clusters does not 
affect our conclusions. The other seven clusters are ro - 
bustly spectroscopically confirmed (jGilbank et alJ f2009) . 

For three of the X-ray clusters (all except WARP 
J1415.1+3612). lMei etafl (|2009D place constraints on z t 
with single-burst models, using a combination of HST 
data, manual morphological measurements, and spectro- 
scopic confirmation where available. They generally find 
formation epochs of zi ~ 2.2, only slightly higher than 
our own measurements. As part of the validation process 
of our pipeline, we compare our measured CMR param- 
eters in these 3 X-ray clusters with previous work, and 
find that they agree within errors. 

For the 4 optical clusters in this study that have ve- 
locity dispersion estimates, we find the distribution of 
velocity dispersions is similar to the X-ray sample, with 
the RCS clusters having 670 < a v < 1080 km s" 1 , and 
the X-ray sample 650 < a v < 747 km s _1 . The number of 

9 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele- 
scope Institute. STSci is operated by the Association of Univer- 
sities for Research in Astronomy, Inc. under the NASA contract 
NAS 5-26555. The observations are associated with program 9425. 
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Fig. 2. — Radial Dependence of E/SO Fraction. X-ray selected 
clusters (solid) show a larger fraction of color-selected E/S0 galax- 
ies than optically-selected (dotted) clusters, and both fractions ap- 
proximately converge to the GOODS fields values (gray bar) by 
0.5 h,- 1 Mpc. Best-fit regression lines are overplotted as well, and 
the intercept and slope (a, 0) are shown in the legend. The X-ray 
fraction at r < 0.4fe — 1 Mpc is consistent with that measured in 
UPostman et al.ll2005t) . 

E/SO galaxies (Table 1) suggests that the clusters in the 
X-ray sample are as rich as the richest optical clusters, 
and a f-test reveals no significant difference between rich- 
nesses of the two samples. Additionally, the fact that we 
have run the same automated pipeline on all the clusters 
rules out systematic differences that might be expected 
in measurements conducted in different studies, or in dif- 
ferent bandpasses. 

The observed difference in red-sequence scatters be- 
tween the two samples therefore suggests that the galaxy 
populations in the clusters themselves are different. To 
test this further, we measure the composite red-sequence 
bulge-dominated (E/SO) fraction in each sample as a 
function of cluster-centric radius, shown in Figure 2. 
Both samples reach the GOODS field value by 0.5/i _1 
Mpc. Overplotted are best-fit regression lines estimated 
using linmix_err (Kelly 2007), which both exhibit a sig- 
nificant correlation of the fraction with radius. The re- 
gression parameters themselves (Figure 2, legend) dis- 
agree at the ~ 2cr level. This suggests that E/SO galax- 
ies form a larger fraction of the cores of X-ray clusters, 



and fall off faster than the optical clusters before reach- 
ing the field value. At lower redshifts the development 
of the morphology-density relation, particularly through 
the build-up of SO galaxies on the red sequence, appears 
correlated with cluster age ([Dressier et alJ[l997l ) . and so 
these differences also suggest that the optical clusters are 
younger. We note however that the radii in Figure 2 are 
metric, and it is not clear that these trends would be 
mitigated by the use of an aperture scaled to the virial 
radius. The similarity of the X-ray and optical velocity 
dispersions is, however, evidence that the virial radii are 
similar. 

The measurements presented in this study have conse- 
quences for studies of cluster evolution and galaxy for- 
mation. Distant optically-selected clusters likely repre- 
sent a younger state of massive halos compared to their 
X-ray selected counterparts. If cluster galaxies from the 
most highly-biased regions of the Universe formed earlier 
in the X-ray sample, as suggested by our analysis, one 
might be tempted to conclude that the clusters them- 
selves collapsed earlier, and that the presence of a hot 
ICM reflects this more evolved state. In this scenario, X- 
ray selected clusters only represent a subset of the most 
massive halos at z ~ 1, and optical samples fill out the 
part of the mass function that is younger, and not X-ray 
bright. However, absent X-ray data for the full optical 
sample, it is difficult to draw any definite conclusions on 
the relationship of cluster dynamical evolution to X-ray 
luminosity and the red-sequence. 
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TABLE 1 

Clusters and CMR Parameters Used in this Study 



Name 1 


N E/S0 




(i — z) slope 


(i — z) scatter 


^0.5 2 


RCS J2319. 8+0038 


39 


0.91 


-0.066 + 0.012 


0.054 ± 0.006 


57.6 ± 1.3 


RCS J0439.6-2904 


23 


0.95 


-0.053 + 0.019 


0.042 ± 0.008 


23.3 + 0.9 


RCS J1511.0+0903 


12 


0.97 


-0.035 + 0.053 


0.055 + 0.020 


9.4 + 1.0 


RCS J0221.4-0321 


18 


1.02 


-0.042 ± 0.073 


0.056 ± 0.025 


26.9 ± 1.3 


RCS J0220.9-0333 


15 


1.03 


-0.010 + 0.022 


0.046 + 0.012 


14.9 ± 1.4 


WARP J1415.1+3612 3 


12 


1.03 


-0.018 + 0.020 


0.022 ± 0.009 


30.8 + 0.9 


RCS J2345. 4-3632 4 


16 


1.04 


-0.050 ± 0.000 


0.051+0.011 


24.5 ± 1.5 


RCS J2156. 7-0448 4 


7 


1.07 


-0.050 ± 0.000 


0.098 + 0.041 


13.6 ± 1.4 


RCS J0337.8-2844 5 


11 


1.10 


-0.108 + 0.090 


0.095 + 0.032 


18.6 ± 1.3 


RDCS 0910+5422 6 


27 


1.11 


-0.044 ± 0.013 


0.044 ± 0.009 


21.5 + 0.8 


XLSS J0223.0+0436 


18 


1.22 


-0.043 ± 0.022 


0.048 ± 0.016 


13.0 + 1.1 


RDCS J1252. 9-2927 6 


19 


1.24 


-0.038 + 0.011 


0.027 + 0.005 


21.5 + 0.7 


RDCS J0848. 9+4452 6 


23 


1.27 


-0.041 + 0.009 


0.034 ± 0.008 


19.4 + 0.5 



Number of E/S0 galaxies (N^/go), cluster spectroscopic redshifts (z S pec), observed 1775 — zgso CMR parameters. 

2 -/V0.5: The GOODS background-subtracted number of galaxies within 2a of the best-fit line (a is our reported scatter) brighter than 
0.3L* at the cluster redshift, and withi n O.Sfe -1 Mpc. 

3 WARPS1415 X-ray detected cluster ijPerlman et a,l.ll200^ l. 

4 A fixed slope is used to measure the scatter, the best-fit slope is sensitive to the color cuts used. See text. 

5 Spectroscopic confirmation pending. 

6 See Mei et al., 2009 



